Delays, connection topology, and synchronization of coupled chaotic maps 
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We consider networks of coupled maps where the connections between units involve time delays. 
We show that, similar to the undelayed case, the synchronization of the network depends on the 
connection topology, characterized by the spectrum of the graph Laplacian. Consequently, scale- 
free and random networks are capable of synchronizing despite the delayed flow of information, 
whereas regular networks with nearest-neighbor connections and their small-world variants generally 
exhibit poor synchronization. On the other hand, connection delays can actually be conducive to 
synchronization, so that it is possible for the delayed system to synchronize where the undelayed 
system does not. Furthermore, the delays determine the synchronized dynamics, leading to the 
emergence of a wide range of new collective behavior which the individual units are incapable of 
producing in isolation. (http://link.aps.org/abstract/PRL/v92/el44101 ) 
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Recent years have witnessed a growing interest in the 
dynamics of interacting units. Particularly, a large num- 
ber of studies have been devoted to synchronization in 
a variety of systems (see lj and the references therein) , 
including the coupled map lattices introduced by Kaneko 

0. Usually, such systems have been investigated under 
the assumption of a certain regularity in the connection 
topology, where units are coupled to their nearest neigh- 
bors or to all other units. Lately, more general networks 
with random, small-world, scale- free, and hierarchical ar- 
chitectures have been emphasized as appropriate models 
of interaction 0, 0, S S • On the other hand, realistic 
modeling of many large networks with non-local inter- 
action inevitably requires connection delays to be taken 
into account, since they naturally arise as a consequence 
of finite information transmission and processing speeds 
among the units. Some numerical studies have regarded 
synchronization under delays for special cases such as 
globally coupled logistic maps || or carefully chosen de- 
lays In this Letter we consider synchronization of 
coupled chaotic maps for general network architectures 
and connection delays. Because of the presence of the 
delays, the constituent units are unaware of the present 
state of the others; so it is not evident a priori that 
such a collection of chaotic units can operate in unison, 

1. e. synchronize. Based on analytical calculations, we 
show that this is indeed possible, and in fact may be fa- 
cilitated by the presence of delays. Moreover, while the 
connection topology is important for synchronization, the 
delays have a crucial role in determining the resulting col- 
lective dynamics. As a result, the synchronized system 
can exhibit a plethora of new behavior in the presence of 
delays. We illustrate the results by numerical simulation 
of large networks of logistic maps. 

We consider a finite connected graph T with nodes 



(vertices) i, writing i ~ j when i and j are neighbors, 
that is, connected by an edge, and with the number of 
neighbors of i denoted by rii. On T, we have a dynamical 
system with discrete time t £ Z, with the state Xi of i 
evolving according to 



/ \ 



xt(t+l) = f(Xi(t))- 



r))-f(xi(t)) 



V iU ) 

(1) 

Here / is a differentiable function mapping some finite 
interval, say [0,1], to itself, e £ [0,1] is the coupling 
strength, and r € Z + is the transmission delay between 
vertices. A synchronized solution is one where the states 
of all vertices are identical, 

Xi{t) = x{t) for alii. (2) 

Thus x(t) satisfies 

x(t + 1) = (1 - e)f{x(t)) + ef(x(t - r)). (3) 

In order to investigate the stability of the synchro- 
nized solution (see |l| for a description of the stabil- 
ity concept employed here), we consider orthonormal 
eigenmodes u k of the graph Laplacian Ar, defined by 



(A r «)< 



Vi), with corresponding 



eigenvalues — Afc. The spectrum of Ar reflects the un- 
derlying connection topology, with the zero ei genv alue Ao 
corresponding to the constant eigenfunction [13j. Since 
the eigenf unctions yield an L 2 -basis for functions on T, 
it suffices to consider perturbations of x of the form 



Xi(t) = x(t) + 5a k (t)uf(t) 



(4) 



for some small 5, and for k > 0, that is non-constant ones. 
The solution x is stable against such a perturbation when 
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&k{t) —> for t — > oo. Expanding about 6 = yields 

a k (t + l) = (l-e)f'(x(t))a k (t) 

+ e(l- X k )f'{x(t-T))a k {t~T). (5) 



The sufficient local stability condition is 

. lim ^™<°- 

T^ooT \a k {0)\ 



(6) 



In the case without delay, that is, t = 0, this is rewritten 



T-l 



<**(* + 1)1 



'-» r "°££ !«*(*)! 



< o 



(7) 



T-l 



and, by (0, becomes 

logll-eAfel + ^i^logl/'^t))! <0, (8) 
that is 



t=o 



|e"(l-eA fc )| < 1 



(9) 



where fi = liniT^oo ^ X)t Jo 1°& i s tne Ly a ~ 

punov exponent of /. 

Essentially the same reasoning works for the case of 
non-zero delay r. Depending on whether |afc(t)| is larger 
than \a k (t— r)| or not, in J2J, we either keep the quotient 
\a k (t + l)/a k (t)\ or replace it by \a k (t + l)/a k (t - r)|. 
This reduces the number of factors in J7J) by r, but since 
we also get correspondingly fewer factors in the term 
leading to the first term in (JSJ, the left hand side of (JSJ) 
is changed by a positive multiplicative factor which does 
not affect the inequality. The terms get slightly more 
complicated, as instead of 



(l-e\ k )f'(x(t))a k (t) 
yielding the two terms in ijSJ), we now have 



(10) 



{l-e)f(x(t))a k (t)+e(l-\ k )f(x(t-T))a k (t-T) (11) 

from JSJ. When comparing l|l(J|l and (|llfl . we see that in 
(|llfl . we may get additional partial cancellations due to 
different signs of f'(x(t)) and f'(x(t — r)) or a k (t) and 
dk {^ — t), so that the absolute value can be significantly 
smaller, thereby making the synchronization condition 
easier to achieve. On the other hand, in a chaotic regime 
the values of f'(x(t)) and f'(x(t — r)) should be essen- 
tially uncorrelated so that we still obtain the Lyapunov 
exponent [i as the relevant parameter coming from /. A 
detailed analysis will appear elsewhere, but the foregoing 
ideas suggest that synchronization in the delayed case is 
at least not more difficult, but potentially easier (due to 
cancellations in dJ) to achieve than in the non-delayed 
one. Moreover, the connection topology, as reflected by 
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FIG. 1: Synchronization of coupled logistic maps for differ- 
ent values of coupling strength e and connection delays r, for 
the cases of (a) scale-free, (b) random, (c) small-world, and 
(d) nearest-neighbor coupling. The grayscale encodes the de- 
gree of synchronization, with black regions corresponding to 
complete synchronization. 
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the eigenvalues still plays an important role in syn- 
chronization. 

As an application, we take / to be the logistic map 

/(a;) = px(l - x) (12) 

which possesses a rich dynamical structure depending on 
the value of the parameter p ^(J . Figure ^ shows the 
synchronization regions in the parameter space for sev- 
eral common network architectures at the value p = 4 for 
which / is fully chaotic. The gray-scale encoding repre- 
sents the degree of synchronization of the network after 
an initial transient of 10,000 time steps, as measured by 
the fluctuations <r 2 (t) = YliLi( x i(t) ~ ^W) 2 j where N 
is the size of the network and x(t) = (1/N) YliLi x iif)- 
Thus, <J 2 (t) — ► as t — > oo if the system synchronizes. 
Darker colors in the figure correspond to smaller values 
of a 2 , with black indicating that a 2 (t) < 10~ 25 after the 
transients. The networks used in the simulations have 
the same size, N = 10, 000, and the same number of av- 
erage connections, even though the architectures may be 
different 0]. The effects of the network topology are 
clearly seen in Figure^ Scale-free and random networks 
can synchronize for a large range of parameters whereas 
more regular networks with nearest-neighbor and small- 
world type coupling do not. In this respect, the simi- 
larities between scale-free and random networks (Figure 
^i-b) are noteworthy. A closer inspection reveals some 
common features for synchronizing networks. For strong 
coupling (roughly for e > 0.6) synchronization is achieved 
regardless of the actual value of the delay, as long as 
it is positive. For intermediate coupling in the range 

0. 4 < e < 0.6, the value of the delay becomes decisive for 
synchronization. In this range one may also observe long 
transients, clustering, and on-off intermittency pj, in the 
gray regions where e is slightly below the synchronization 
value. There are also smaller regions of synchronization 
that exist for weaker coupling (0.15 < e < 0.20) and only 
for odd delays. Note that for zero delay synchronization 
can occur only for a rather limited range (e > 0.85). Also, 
the small region of synchronization in Figure Q>d occurs 
for nonzero delay. Therefore, the presence of delays can 
indeed facilitate synchronization. 

While the connection topology is important for the 
synchronizability of the coupled system, connection de- 
lays are significant in determining the resulting synchro- 
nized dynamics. In particular, there is an important dif- 
ference with the undelayed case, for which 10) reduces to 

x(t + l)=f(x(t)); (13) 

1. e., when the synchronized system can exhibit nothing 
different from the exact dynamics of the individual iso- 
lated unit. Hence, an important implication of connec- 
tion delays is the possibility of the emergence of new col- 
lective phenomena. Indeed, the solutions of J3J exhibit a 
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FIG. 2: The bifurcation diagram of the synchronized solution 
of coupled logistic maps for several values of r. At each value 
of p, 200 iterates of (J^J are plotted after inital transients. The 
value of e is 0.8. 



much richer range of dynamics when r is nonzero. Fig- 
ureyjshows the bifurcation diagram for several values of r 
when / is given by l|12jl. Of course, for r = the familiar 
bifurcation diagram of the logistic map is obtained, dis- 
playing the period-doubling route to chaos. By contrast, 
when r > Neimark-Sacker type bifurcations are 
prevalent, immediately resulting in high-period solutions 
followed by more complex behavior. Other values of odd 
and even delays result in pictures similar to the cases 
t = 1 and r = 2, respectively. Since Neimark-Sacker 
bifurcations cannot arise in one-dimensional maps, the 
difference with the undelayed case is fundamental. Sim- 
ilarly, a consideration of the Lyapunov exponents shows 
that the character of the chaotic attractor is also dif- 
ferent. For instance, as seen in Figure the synchro- 
nized system can have two positive Lyapunov exponents 
when connection delays are present, thus exhibiting hy- 
perchaos. 

Connection delays make the coupling strength e an im- 
portant factor in the dynamics of the synchronized solu- 
tions. To demonstrate, we take the logistic map <|12[) 
with p — 4, for which it has a chaotic attractor. Thus, 
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FIG. 3: The Lyapunov exponents of the synchronized solution 
for r = 1 and e = 0.8. 
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FIG. 4: The synchronized solution at different coupling 
strengths for r = 1 and p = 4. At each value of e, 300 it- 
erates of © are plotted after inital transients. 

© is chaotic for all e when r = 0. By contrast, when 
t = 1 the attracting solutions of display a wider vari- 
ety depending on the value of e, as depicted in Figure 01 
Chaos is interrupted by windows of stable periodic solu- 
tions, for instance near e = 0.68 and s = 0.71. There is 
a large window of period-5 solutions for 0.82 < e < 0.89, 
which also contain other stable solutions of higher peri- 
ods for certain e in this range. A periodic attractor of 
period 3 exists near e = 0.945, which transforms to in- 
termittent behavior at e — 0.95, shown in Figure [SJ All 
these dynamics are exhibited by the coupled system QJ, 
since Figure ^ implies that for these parameter values 
the system is synchronized (provided it has an appro- 
priate connection topology), so the collective behavior is 
described by J3J ■ Clearly, delays can make the dynamics 
of the synchronized system quite sensitive to the coupling 
strength, a feature that is absent in undelayed networks. 

Networks with connection delays arise naturally in 
many areas of science, including biology, population dy- 
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FIG. 5: Intermittent behavior of synchronized logistic maps 
for p — 4, e = 0.95 and r = 1. 



namics, neuroscience, economics, and so on. In neural 
systems, for instance, delays result from e. g. finite ax- 
onal transmission speeds [l2j. Such networks lack an 
intrinsic notion of simultaneity since the present state of 
the system is inaccessible to the constituent units. These 
networks are nevertheless capable of operating in syn- 
chrony, even when complex units and significantly large 
delays are involved, as our findings indicate. In fact, net- 
works with delays can actually synchronize more easily. 
This synchronizability property is especially relevant for 
neural networks, circumventing the difficulties in estab- 
lishing a concept of collective or simultaneous informa- 
tion processing in the presence of delayed information 
transmission [15j. Furthermore, the delays shape the dy- 
namics of the synchronized system, leading to the emer- 
gence of a variety of new dynamics which the individual 
units are not capable of producing . The observation 
that a wide range of different behavior can be accessed by 
varying the coupling strength has important implications 
for neural networks, where synaptic coupling strengths 
can change through learning. This interesting connec- 
tion provides additional motivation for investigating the 
role of delays in complex adaptive systems. 
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